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Abstract 

A novel centrifuge set-up for the study of unsaturated flow charac- 
teristics in porous media is examined. In this set-up, simple boundary 
conditions can be used, but a free moving boundary between unsaturated- 
saturated flow arises. A precise and numerically efficient approximation is 
presented for the mathematical model based on Richards' nonlinear and 
degenerate equation expressed in terms of effective saturation using the 
Van Genuchten-Mualem ansatz for the soil parameters in the unsaturated 
zone. Sensitivity of the measurable quantities (rotational moment, center 
of gravity and time period to achieve quasi steady state) on the soil pa- 
rameters is investigated in several numerical experiments. They show that 
the set-up is suitable for the determination of the soil parameters via the 
solution of an inverse problem in an iterative way, excluding the saturated 
hydraulic conductivity. For this parameter, an existing simple centrifuge 
set-up is repeated and augmented with transient measurements. 
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1 Introduction 



To predict the flow and solute transport in soils one needs the soil hydraulic 
properties in terms of soil parameters. Once determined, these parameters can 
be used as input data in the governing mathematical modeL For unsaturated 
flow, this model is given in terms of the saturation and the pressure head in 
Richards' equation (see below), which is a nonlinear and degenerate parabolic 
equation. Furthermore, when part of the sample is saturated, free boundaries 
between the saturated zone and the partially saturated zone arise, as well as 
between the dry and the partially saturated zone. This is a major problem 
for many modeling approaches, leading to experimental set-ups that avoid the 
formation of these boundaries. 

The soil retention and hydraulic permeability functions linking the satu- 
ration and pressure head for unsaturated flow are expressed using the Van 
Genuchten-Mualem ansatz by means of soil parameters. Measuring these soil 
parameters is usually time consuming and tedious, especially for low conductive 
porous media. Several set-ups based on centrifugation have been proposed to 
obtain a large acceleration of the processes involved, see [2 |4j [5j [3 13] and cita- 
tions therein. These techniques have several disadvantages. An approach which 
aims for a steady-state flow regime inside the centrifuge, [US], requires expen- 
sive and/or complex apparatus, and obtains only a few water content versus 
conductivity measurements per run. Also transient set-ups based on keeping a 
top boundary at a fixed prescribed setting, |3j, are expensive. The quasi-steady 
centrifuge (QSC) method, [T], is a much simpler technique (a slowly emptying 
reservoir at the top that is refilled when needed), but requires that the cri- 
terion for steadiness of flow through the sample is relaxed, leading to higher 
uncertainty in the obtained results. 

The alternatives for determining conductivity with a steady-state flow, cou- 
ple a transient flow with parameter estimation techniques, see eg. (21 [5]. In 
this way, the conductivity and retention curve can be determined inversely over 
a large saturation domain. These methods require experiments of some state 
variables which relate to the conductivity. One-step or multi-step outflow meth- 
ods are common in column experiments. Pressures or suctions are applied at 
the top or the bottom of the sample, while the outflow, and optionally pressure 
head within the sample, are measured. These measurements are then used to 
estimate the hydraulic parameters. This technique is transferred to the cen- 
trifuge device in [8]. However, measuring precise outflow rates is problematic 
in a centrifuge. Therefore, their method for determination of the soil parame- 
ters is based on transient water content measurements (coming from electrodes 
installed in the specimen). This can be used in an equilibrium analysis as well 
as in a transient analysis. In this setting, the specimen is fully saturated at the 
beginning, and saturation is maintained at the outflow via a short lip at the 
end of the container. Good results are obtained, but there remain some dis- 
advantages to this technique: there are few measurements close to saturation, 
leading to a high error in the prediction of the conductivity close to saturation, 
the sample needs to be disturbed to introduce the electrodes, and there is a 
very long waiting time in order to achieve the equilibrium when the equilibrium 
analysis approach is used. 

The main goal of this manuscript is to develop a precise numerical method 
enabling to determine the soil parameters (via solution of inverse problem) in a 
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very simple way requiring very cheap measurements. 

In this article we consider any partially saturated sample which is sealed 
at the right boundary (from the center of centrifuge) and has no inflow at the 
left boundary. The amount of infiltrated water is fixed and would be important 
initial information, but in the fact, we can drop this requirement and consider 
the initial water content as an additional parameter requiring determination. 
The only measurements required in our setting is the rotational momentum 
and the center of gravity of the sample in a set of equilibria corresponding 
to some predetermined rotational speeds. This can be replaced however by 
with transient measurements of these variables allowing to reduce the centrifuge 
operating time. These measurements are sufficient due to the fact that the 
saturation profile at the equilibria is not depending on the initial distribution 
of water in specimen, but only on its amount, which, when the right boundary 
of the sample is sealed, is identical in all equilibria. 

To use this procedure, we have to face serious difficulties in the numerical 
modeling. The main one is that if the right side of the sample reaches effective 
saturation, an interface between partially saturated zone and saturated zone 
appears. This boundary is very difficult to control numerically, causing problems 
with the mass balance conservation which is very important in this set-up. 

To reach the equilibrium is an infinite asymptotic process, but after some 
time (eg. 1-3 days for low conductive material) the change of the rotational 
momentum and of the center of gravity can no longer be measured. At that 
moment, the rotational speed is increased, and the system moves towards a new 
corresponding equilibrium. Note that even when equilibrium was not reached 
and a small error is present in the measurements of the rotational momentum 
and the center of gravity, this will not influence the error at the higher equilib- 
rium level. This error depends only on the running time of centrifugation at the 
actual rotational speed. The differences between applied rotational speeds are 
chosen in such a way that that the differences in outputs (rotational momentum 
and center of gravity) are technically well distinguishable. 

Next, the soil parameters and eventually the amount of originally infiltrated 
water, can be determined by minimizing a cost functional expressing the dis- 
tance between the measured and the computed output, e.g., with the Levenberg- 
Marquard method. In this manuscript we only show the sensitivity of the out- 
puts on the parameters, which is a good indication of how an inverse method 
will perform. In fact, the correct determination of soil parameters substantially 
depends on the sensitivity and the errors in practical measurements. 

We can also use the transient measurements of the outputs, for example the 
period needed to cross from one equilibrium to another. This would give us 
additional information useful in the determination of the soil parameters. 

The advantage of the above approach is that the full range of saturation 
values are present in the setup, while preventing outflow means equilibrium can 
be obtained faster. However, due to the set-up, it is clear that the water flows 
from the unsaturated zone to the saturated zone, with no flow occurring in the 
saturated zone. Indeed, we notice that the rotational momentum and center of 
gravity in the equilibria approach are not sufficiently sensitive on the "saturated 
hydraulic conductivity" . Therefore, we repeat in Section H] also centrifugation 
with a fully saturated specimen which substantially increases sensitivity on this 
soil parameter, see [5]. In this setup the sample remains saturated and the water 
table head in the top reservoir is measured over time and related to the saturated 
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hydraulic conductivity. It is in other words the same as a column experiment 
where the dropping head at the top is measured instead of the outflow, and the 
head at the bottom is fixed instead of the head at the top. 

In our numerical method we reduce the mathematical model to a system 
of ODE using method of lines (MOL), which has already been successfully ap- 
plied to Richards' equation in e.g. Our main contribution is in correctly 
handling the moving free boundary. The obtained system can be solved with 
ODE solvers for stiff systems and stiff DAE systems (system of ODE and al- 
gebraic equations). The numerical experiments demonstrate that this solution 
method is numerically very efficient (precise and economical) . Note that reach- 
ing equilibrium for a high rotational speed can require a week or even a month 
of centrifugation for low conductive materials, whereas the computational time 
is only a few seconds. Therefore this method is a good candidate for solving 
the inverse problem, as this requires many computations of the forward prob- 
lem. The numerical method can be successfully applied in other centrifugation 
settings (concerning control of the inflow, or control of the outflow) as, e.g., in 



The research performed in this article was done for a preliminary study 
for the feasibility of a new centrifuge device for low conductive materials that 
should be cheap to operate. The results presented here indicate that an ap- 
proach where only rotational moment and gravity center is needed, is workable. 
This would remove the need to add experimental apparatus to the centrifuge 
measuring pressure heads, water contents or outflow rates during the operation. 
The authors intend to investigate a different set-up in a future article, in order 
to determine an optimal candidate, as well as collaborate to create a prototype 
centrifuge. 

In Section [2] we present the mathematical model, giving specific attention to 
the movement of the free boundary. In Section |3] the numerical method based 
on the MOL approach is given, while in Section U the approach to determine 
the saturated hydraulic conductivity is explained. We finish in Section [5] with 
several numerical experiments showing the sensitivity of the output parameters 
on the soil parameters. 

2 Mathematical model 

We consider a one dimensional model for a partially saturated sample in the 
form of a tube. The tube starts (top or left boundary) at the distance r = 
ro from the center of the centrifuge and ends at the distance r — tq + L. 
The right boundary of the specimen is isolated. Flow in porous media under 
centrifugation is modeled by Darcy's equation in the saturated region and by 
Richards' equation in the unsaturated region (see, e.g., [H],[3])- So 



00- 




(1) 



in the saturated region, and 
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in the unsaturated region. Here, h is the piesometric head, 9 the saturation of 
the porous medium, w the angular speed of rotation (in radians per second) , Kg 
the hydrauhc conductivity in the saturated region, g the gravitational constant 
and the function k{9) describes the hydraulic conductivity in the unsaturated 
region. Denote by m = j—^ the effective saturation, where 9s is the volumetric 
water content at saturation and 9r is the residual volumetric water content. We 
have u G (0, 1), since 9 e {9s, 9r)- The soil hydraulic properties are represented 
by empirical expressions (see [7]) 

^ /ie(-cx),0) (3) 



(1 + (7/1)")" ' 

k{u) = Ksu^'^[l - (1 - (4) 

where m = I — 1/n, n > 1 and 7 — —(2^/™ — l)^^'"//if, are empirical soil 
parameters, hb is called the bubbling pressure . 

It is possible to rewrite the flow in unsaturated form as 



dtu — dr D(u)drU k{u)r ] (5) 

.9 



where 



[n - 1)7(6/^ - 9r) 



X [1 - (1 - (6) 



Equation ([5]) is strongly nonlinear and degenerate. We note that D{Q) = 
0, -D(l) = 00. Equilibria at the high rotational speed can be expected to have 
a fully saturated zone (supposing the initial amount of infiltrated water is suf- 
ficiently large), which appears at the right boundary and of which the front 
evolves to the left of the specimen (under non-decreasing rotational speed) . We 
denote the position of this interface by s{t). This saturated zone is governed by 
Darcy's equation, but s[t) is unknown and time dependent. The time evolution 
of s(t) is difficult to compute. The dynamics of this region is linked with the 
(finite) interface flux qi 



qi = — \ D{u)drU k{u)r 

V 9 / r=s{t) 

and based on a mass balance argument we can expect 

sit) = -qt- 

Unfortunately, we cannot use this model for the determination of the time evo- 
lution of s{t), since at r = s{t) it holds u = 1 and £'(1) = oo. Consequently, 

drU\r=s{t) = 0. 

If we transform Richards' equation in terms of the head, we obtain 



ds{h)dth = k^dr 



^2 



k(h)drh k(h)r 

9 



(7) 
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Figure 1: k{h) and 100 x for n = 2.81, 7 = -0.0189 



with fco = -g^zfg~: where kok{h) is the hydrauhc conductivity function, 

(7/1)"-! 



k{h) = 



1 



1 



(1 + (7/i)")™/2 (1 + (7/1)")'- 

and the specific moisture capacity function ds{h) — du/dh is given by 

(7/1)"-! 



4(M = -7(n-i) 



(1 + (7/1)")!+'' 



We can see that k{h) — s> 1 for h — > 0. In Fig. [T] we present the graph of the 
functions k{h) and 100 ds{h) for G (—200,0), and parameter vahies Kg = 
2.4 10~^, n = 2.81, 7 = —0.0189. As we can see also equation (O degenerates 
at h = 0. This has to be taken into account when saturation becomes 1 at 
the right boundary of specimen. After this moment, t ^ ti, the mathematical 
model must be changed to reflect the physical phenomenon. At the right hand 
side of the (isolated) specimen appears a saturated zone with an interface s{t) 
moving from the right boundary to the left. The flux at the interface s{t) is 
equal to —s(t), but also in this pressure-head form of Richards' equation it is 
difficult to approximate correctly dth\x=s{t), which leads to a significant error 
in the mass balance. 

Therefore to determine the interface s{t), we will consider the algebraic 
equation 



ro+L 



u{h{x,t)) dt + L - s{t) ^ AU, s(0) = L, 



(8) 



where is the amount of infiltrated water (which remains constant during the 
centrifugation). This condition refiects the global mass balance in the specimen 
and does not suffer from a fiux approximation at r = tq + s{t). 

Then mathematical model ([T]) only needs to be solved over the interval r G 
{ro^ro + s{t)) with right boundary condition h{ro + s{t)) = for all t. We 
approximate this mathematical model in the next section. 
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3 Numerical method 



For the output parameters that wiU be measured (gravity center and rotational 
momentum), there is no need to model the head in the saturated zone, as we 
consider the compressibility of water to be negligible. The numerical approxi- 
mation of mathematical model ©-([SI) results in a coupled system of a PDE and 
an algebraic equation. Moreover, the solution domain is a moving region, with 
unknown interface s(t), which has to be determined. 

We shift Q to the domain r G {0,s{t)) and use the fixed domain transfor- 
mation y — -j^. This gives 



d,{h)[dthiy,t)~y^dyh 



kn 



1 



dy [ k{h)dyh - k{h)— (ro + ys{t)) 



(9) 



Consider the space discretization 

yo < < • • • < 2/i < 



< 2/7V = 1, 



and ao = 0, ai yi — j/i-i, i — l,...,iV and integrate ^ over li : 
(2/1-1/2,2/1+1/2) for i = 1, . . . ,iV - 1 where y,-i/2 ■= iUi + y-*-i)/2, y^+i/2 ■ 
{yi + yi+i)/2. 

We denote by hi(t) w h{yi, t), Vi = 1, . . . , A^ — 1, and approximate dth{y, t) 
hi{t) in the interval 7^. We approximate 



'^y^\y=yi+i/2 



ht+ijt) - h,(t) 
ai+i 



and similarly we approximate 9yh\y=y._-^^2 and denote it by : d^hi. Let £(z; yi) 
be the second order Lagrange polynomial crossing the points (j/i-i, /li-i), {yi, hi) 
and (j/i+i, /li+i). We use the abbreviation A:i+i/2 := ^(^j/i+1/2) Then, our ap- 
proximation of ^ (based on finite volume type approximation) at the point 
y = yi reads as follows 



\ s dz 



ko- 



z=yi 
2 1 



i+1/2 



^i-l/2'9 hi 



— {h+i/2{'rQ + sy^+i/2) - fcj-i/2(ro + syi-1/2)) 



(10) 



for i = 1 , . . . , — 1 . We add the corresponding equation at the point yo taking 
into account that the flux is zero there. In a similar way as in (fTO)) (following 
the finite volume type of approximation) we obtain 



2 1 
ds{ho)ho = ko—-T 
all) 



— (fci/2(r-o + syi/2)) 



(11) 
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At the point ^at = 1 we have h^it) — 0, so no additional equation is considered. 
We approximate the amount of water using the trapezoidal rule for the 
integration. Define 



Q{t) uo ai/2 + aAr/2 + ^ 



2 



1 




' (l + (7/i,)")'"' 
Then, system (|10p - pT|) will be completed by the algebraic equation 



= L-s{t)[l-Q{t)]-M^. 



(12) 



This algebraic equation is used instead of an ODE equation that models s{t). 
System (fTUl) - ([T^ is degenerate and is of the form 



where z — (ho,hi, ...,hN^i, s). The last equation of this system is just (|12p. 
This system can be readily solved, e.g., by the solver "odel5s" in MATLAB® 
or the "ida" solver of the Sundials package. 

As is usual with these solvers, some regularization in ([T2|) is needed as well 
as a tuning of the space discretization. Most important is to have a "good" 
starting point. 

If the equilibria have the property < 0, then no interface appears. It is 
then needed to set s{t) L in the previous mathematical model and replace 
algebraic equation (jl2p by an ODE equation for /i^v which will be similar to 
((TT|l . Successively increasing the rotational speed of the centrifuge increases the 
head at the right boundary. The model remains in the state where s(t) := L up 
to the point when h{N) — 0, at which point the computation is automatically 
halted. The full model (fTU)) - ([T^ is used onwards to compute the equilibrium 
states. 

In the numerical equilibrium experiments (see Section [5l ExD. I5.5l and l5.7p it 
is observed, as expected, that the values of the rotational moment Mr and the 
center of gravity Gc are not very sensitive on the important Kg parameter. Also 
the transient experiments where the time sections between different equilibria 
are measured, are not very sensitive. The saturated conductivity Ks can only 
be determined from measurements of Mr, Gc that are accurate up to 3 digits. 
Therefore, another method must be used for the determination of A'.;, which we 
present in the next section. 

4 Centrifugation of saturated sample 

For the determination of the saturated conductivity we propose to use the device 
put forward in [2, with the addition of allowing for transient measurements. We 
specifically use the ability to measure when a reservoir has completely drained 
out, combined with the measurements of the rotational moment. 

We consider a saturated sample placed in the region (ro,ro + L) of the 
centrifuge. A water reservoir is placed in front of the sample, r e (rp — IotTq) 



M{t,z)z{t)^ f{t,z) 



(13) 
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and a head of is maintained at the bottom rQ + L via a Up (alternatively one 
could place an empty reservoir r S (tq+L, rQ+L+d) into which the water flowing 
out from the sample is collected). The water from the front reservoir will flow 
through the sample and drain out to the right. We proceed the centrifugation 
up to the moment Tg which is defined as the moment when the front reservoir 
is empty. By means of this setting the measurements of rotational momentum, 
resp. Te, will be very sensitive on the model parameter Ks- The information of 
only the end time Te might even be sufficient to determine Kg . 

The mathematical model is very simple, since we have a saturated sample 
during the centrifugation - see ([1]), with time dependent Dirichlet boundary 
condition. The flux of water is given by 



qF{t) = -K,l^drh-—rj, (14) 

and is constant along all the sample. The head at the left boundary r = tq 
depends on the water level £{t) (in the reservoir) and is under centrifugation 
given by 

2.9 

The head at the outflow boundary is zero, i.e., h{rQ + L) ^ 0. Hence, by simple 
integration of ([T} in the saturated region we obtain the solution 



h{ro + x,t) = — [£{t){2ro - m) + 



L J 2roL 
where x £ (0,i). Inserting ([15]) in (ITU) we obtain 



(15) 



qpit) ^Ks— [L^ - £{tf + 2ro(L + £{t))\ (16) 

Due to the mass balance argument we obtain the governing ODE for the time 
evolution of the water level: 

with £(0) = /q and l{Te) = 0. Solving this ODE we obtain the relation between 
Tg and Kg, whereas £{t) fully determines the change of the rotational moment 
Mr{t) over time. The numerical experiments for this set-up is presented in Sec. 
[Ql below. 



5 Numerical experiments 

For all the experiments we use as data rg = 10, L = 10, uj = 30, Kg = 2.4 10~^, 
9r — 0.02, 9s = 0.4, 7 — —0.0189, n — 2.81, except where noted otherwise 
or where sequences are compared to investigate the sensitivity of the set-up 
on the parameters. If not mentioned otherwise, a uniformly distributed space 
discretization with iV = 40 grid points is used. 
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Figure 2: Evolution to equilibrium starting from a uniform water distribution, 
w = 30 

5.1 Independence of initial data 

In the first numerical experiment we demonstrate that the equilibrium depends 
only on the rotational speed and the amount of infiltrated water. In Fig. [2] 
we show the time evolution of a set of unsaturated profiles over 10 equidistant 
time sections in the interval (0,500.000s). From this figure it is clear that at 
the end time the profile is very close to the equilibrium (of which the precise 
value is obtained in infinite time). Moreover, a long time before the end of 
the centrifugation at t = 500.000s, the saturation profiles do not change visibly 
anymore. 

The independence of the initial profile is clear in Fig. [3] where the only change 
is in the starting saturation profile (which is now lineary decreasing), keeping 
the amount of infiltrated water identical. The same equilibrium as in Fig. [2] is 
approached. 

5.2 Rotational speed 

We investigate the sensitivity on the rotational speed w. For this we monoton- 
ically increase w in a stepwise way, from a; = 10 up to w = 70 with increment 
Aco = 5, resulting in 13 saturation profiles. At each rotational speed we reach 
the corresponding "good approximated" equilibrium. The centrifugation from 
one equilibrium to the next one takes 800.000s « 9days. In Fig. 2] we present 
the saturation profiles of the equilibria, together with one extra profile corre- 
sponding to the right boundary becoming fully saturated. The corresponding 
heads are presented in Fig. [5j 

The extra profile is the sixth curve, which occurs during w = 35 and has value 
— 0. At that moment, the ODE solver is stopped, the system is changed 
to the system with a moving boundary, Eq. ©-([H]), and the solver is restarted 
with initial condition given by the last potential head (resp. saturation) profile. 

The values at the equilibria of the rotational momentum Mr, the center 
of gravity Gc, and total water content are collected in Table [TJ These 
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Figure 3: Evolution to equilibrium starting from a lineary decreasing water 
distribution, a; = 30 
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Figure 4: Equilibrium saturation profiles for w = 10 + j 5, j — 0,...,12, and 
profile with ft, at = 
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Figure 5: Equilibrium head profiles for lo = 10 + j 5, j = 0,...,12, and profile 
with /ijv — 

are computed from the saturation profiles. Two rotational moments are given, 
Mr.s denotes the starting rotational momentum corresponding to the satura- 
tion distribution obtained from the previous rotational speed, and M^^e denotes 
the rotational momentum at the actual rotational speed in equilibrium. The 
formulas for M^, Gc and at time t are: 



and are all evaluated numerically using the trapezoidal rule. Note that if 
u{t, 1) < 1 then sit) — L. The sensitivity of the measured quantities on the 
changing water content is very good. It remains to determine the contribution 
of the different parameters. 

We note first however that careful analysis of the head profiles in Fig.[S]shows 
that these are not yet parabola at the end times of a run at fixed rotational speed, 
especially toward the left of the sample. This indicates that the equilibrium has 
not yet been reached. We investigate this further in the next Experiment. 

5.3 Reaching equilibrium 

To investigate the head profiles we start in this experiment from the equilibrium 
corresponding to w = 40 and use a rotational speed w = 50. The centrifuge 
normally operates up to Tg = 1.540.000 seconds. At that time, equilibrium for 
a; = 50 is almost reached. We compare 13 values, the starting value, 9 increasing 
time steps (with Atj = ij+i — tj = 2000 2^, j — 1, ...,9), the sensible end time 
step Te — 770 X 2.000s, and 2 extra time steps to investigate the very long 



Mr = 



/•I 1 

/ (ro + sit)z)^u{t, z)dz + -{L^ - s{tf), 
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Table 1: Rotational momentum, center of gravity, water amount for Exp. [521 





Mr,s lO"'' 


Mr,e 10"" 


Gc 


M„ 


10 


0.0467 


0.0481 


4.9999 


3.9999 


15 


0.1081 


0.1126 


5.4979 


4.0020 


20 


0.2002 


0.2118 


5.8722 


4.0041 


25 


0.3309 


0.3498 


6.3012 


4.0057 


30 


0.5037 


0.5288 


6.7097 


4.0087 


35 


0.7037 


0.72307 


6.8220 


4.0110 


40 


0.9307 


0.9747 


7.2727 


4.0141 


45 


1.2338 


1.2534 


7.4323 


4.0137 


50 


1.5474 


1.5649 


7.5494 


4.0133 


55 


1.8936 


1.9089 


7.6349 


4.0130 


60 


2.2718 


2.2853 


7.6984 


4.0126 


65 


2.6821 


2.6941 


7.7469 


4.0123 


70 


3.1245 


3.1352 


7.77845 


4.0120 



Table 2: Rotational momentum, center of gravity, water amount for Exp. 15.31 



time 


M,. e 10-^ 


Gc 







1.5201 


7.2512 


4.0141 


1 


1.5248 


7.2813 


4.0128 


3 


1.5299 


7.3119 


4.0131 


7 


1.5345 


7.3413 


4.0133 


15 


1.5389 


7.3697 


4.0134 


31 


1.5430 


7.3972 


4.0134 


63 


1.5469 


7.4234 


4.0135 


127 


1.5505 


7.4478 


4.0135 


255 


1.5537 


7.4699 


4.0136 


511 


1.5565 


7.4893 


4.0136 


770 


1.5588 


7.5056 


4.0136 


1800 


1.5645 


7.5462 


4.0133 


2300 


1.5649 


7.5494 


4.0133 



time behavior. Fig. [S] shows the obtained saturation profiles, the corresponding 
evolution of the potential head is presented in Fig. [T] As we can see the 
head profile at Tg is still not a parabola in the left side of the profile (low head 
values). The reason for this is that the hydraulic permeability at low head is 
negligibly small, so it takes a very long time to reach the equilibrium. If the 
centrifugation is continued, also this part obtains the required parabolic shape 
associated with the equilibrium. Note however, that the other part of the head 
profile (for higher head values) is changing insignificantly. Therefore, we arrive 
at the conclusion that it makes sense to increase the rotational speed and not 
wait for these lower head values to stabilize. More so in light of Exp. 15.11 

The measured values for the rotational momentum, gravity center, and water 
amount, are given in Tabled The small change between the last two values in 
Table [2] demonstrates that equilibrium is eventually reached. 
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Figure 6: Evolution of the saturation starting from equilibrium at w = 40, with 
w = 50 




space 

Figure 7: Evolution of potential head corresponding with Fig. [6] 
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Figure 8: Equilibrium saturation profiles at w — 20 for n — 1.51; 1.81; 2.11; 
2.41; 2.71; 2.81; 3.01 and 3.31. 



Table 3: Rotational momentum, center of gravity and water amount for Exp. 
l Ol , , 



n 


M^e.lO-" 


Gc 


M„ 


1.51 


0.1887 


5.0736 


4.0043 


1.81 


0.1927 


5.2391 


4.0052 


2.11 


0.2020 


5.6188 


4.0043 


2.41 


0.2068 


5.8096 


4.0062 


2.71 


0.2083 


5.8701 


4.0054 


2.81 


0.2112 


5.9870 


4.052 


3.01 


0.2153 


6.1524 


4.0060 


3.31 


1.5505 


7.4478 


4.0135 



5.4 Dependence on n 

In this experiment we demonstrate the sensitivity of Mr and Gc on the model 
parameter n. We start with a constant saturation u = 0.4 and apply now the 
rotational speed lu = 20. The centrifuge is operated for 800.000s. In Fig. [5]the 
obtained equilibrium profiles are depicted for successively n = 1.51; 1.81; 2.11; 
2.41; 2.71; 2.81; 3.01 and 3.31. In Fig. [3] the corresponding heads are given. 

The resulting values for Mr , Gc and M^ are given in Table [31 and indicate 
a good sensitivity. 

5.5 Saturated conductivity 

We now proceed with investigating the sensitivity of Mr and Gc on the saturated 
conductivity Ks- We consider = h 10"^ for = 0.4; 0.8; 1.6; 2.0; 2.4; 2.8 
and 3.2, start again from a uniform saturation, and apply this time a rotational 
speed of w = 20. In Table H] we present the resulting values of Mr, Gc and M^ 
at the obtained equilibrium. The sensitivity of Mr and Gc on ks is negligible at 
the relative low pressure head that is generated for a; = 20. 
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Figure 9: Equilibrium head profiles corresponding with Fig. |8] 



Table 4: Rotational momentum, center of gravity and water amount for Exp. 
531 = 20 



kg 


M^^e.lQ-'^ 




M„ 


0.4 


0.2049 


5.7286 


4.0052 


0.8 


0.2077 


5.8401 


4.0078 


1.2 


0.2081 


5.8633 


4.0059 


1.6 


0.2087 


5.8684 


4.0098 


2.0 


0.2082 


5.8697 


4.0074 


2.4 


0.2083 


5.8700 


4.054 


2.8 


0.2083 


5.8703 


4.0062 


3.2 


0.2087 


5.8702 


4.0095 
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Table 5: Rotational momentum, center of gravity and water amount for Exp. 
[531 = 50 





Mr.e 10-^ 


Gc 




0.0001 


0.1867 


5.0118 


4.0006 


0.001 


0.1888 


5.0855 


4.0025 


0.01 


0.1998 


5.5221 


4.0059 


0.1 


0.2085 


5.8697 


4.0063 


1.0 


0.2084 


5.8699 


4.0074 


10 


0.2083 


5.8701 


4.0064 


100 


0.2083 


5.8703 


4.0022 



0.9 
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Figure 10: evolution of saturations at equilibrium for uj = 50, kg = 
0.0001; 0.001; 0.01; 0.1; 1; 10; 100 



In Table [5] we present analogous results for the equilibrium obtained at ro- 
tational speed u) = 50, starting from the equilibrium at oj ~ 40, and considering 
a much larger spectrum for ks . To reach the equilibrium from which Table [S] is 
extracted, a centrifugation time of 15 lO^s is used for the values fcs < 1, whereas 
for values fc^ > 1 we have used an operating time of lO^s. 

In Fig. [TU] the equilibrium profiles of the saturation is given corresponding 
with the results from Table [SJ while the corresponding heads are given in Fig. 
im From Fig. [TT] one can see that for the lower values of the head a parabolic 
shape of the profile was still not reached if fcs < 1. Overall, the results clearly 
indicate that another technique must be used to determine K^- 

5.6 Dependence on 7 

We now investigate the sensitivity of Mr and Gc on the last model parameter: 7. 
We again use a rotational speed of a; = 50, starting from the equilibrium position 
dXoj — 35. As values for 7 we consider 7 — —70 10^ with 70 £ (1.59; 2.19) where 
increments of size 0.1 are used. The values of Mr and Gc are listed in Table 
[6] and Table [71 respectively. The corresponding saturation and head profiles 
at time section t = 10^ are in Fig. [TU and Fig. [T31 respectively. In Tables 
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Figure 11: evolution of heads at equilibrium for uj 
0.0001; 0.001; 0.01; 0.1; 1; 10; 100 



50, ks 



Table 6: Rotational momentum Mr 10 ^ for Exp. 



time\7o 


1.59 


1.69 


1.79 


1.89 


1.99 


2.09 


2.19 


1000 


1.5189 


1.5093 


1.5013 


1.4949 


1.4896 


1.4854 


1.4819 


3000 


1.5352 


1.5213 


1.5097 


1.5000 


1.4919 


1.4852 


1.4797 


5000 


1.5438 


1.5278 


1.5144 


1.5030 


1.4935 


1.4855 


1.4788 


10* 


1.5565 


1.5376 


1.5216 


1.5079 


1.4963 


1.4864 


1.4780 


5.10* 


1.5901 


1.5645 


1.5423 


1.5231 


1.5063 


1.4917 


1.4791 


10^ 


1.6058 


1.5777 


1.5530 


1.5313 


1.5124 


1.4958 


1.4812 



Table 7: Center of gravity for Exp. 



time\7o 


1.59 


1.69 


1.79 


1.89 


1.99 


2.09 


2.19 


1000 


7.1117 


7.1031 


7.0948 


7.0872 


7.0805 


7.0744 


7.0693 


3000 


7.1487 


7.1379 


7.1271 


7.1170 


7.1076 


7.0990 


7.0912 


5000 


7.1697 


7.1581 


7.1464 


7.1351 


7.1245 


7.1146 


7.1056 


10* 


7.2026 


7.1903 


7.1775 


7.1649 


7.1527 


7.1411 


7.1303 


5.10* 


7.3012 


7.2898 


7.2763 


7.2617 


7.2466 


7.2314 


7.2167 


10^ 


7.3512 


7.3429 


7.3309 


7.3167 


7.3012 


7.2851 


7.2689 
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Figure 12: Saturation profiles at equilibrium for uj = 50, 70 — 
1.59; 1.69; 1.79; 1.89; 1.99; 2.09; 2.19 



Table 8: Rotational momentum vs. under changing uj at the time level 
= 1000 + j 10^ 





fcl 


fc2 


fc3 


/C4 


fcs 


fce 


fc7 


fcs 


fcg 


to 


0.9526 


0.9531 


0.9536 


0.9538 


0.9540 


0.9555 


0.9564 


0.9570 


0.9576 


h 


1.2123 


1.2169 


1.2196 


1.2213 


1.2228 


1.2305 


1.2339 


1.2358 


1.2370 




1.5049 


1.5147 


1.5200 


1.5237 


1.5265 


1.5401 


1.5450 


1.5475 


1.5490 




1.8301 


1.8458 


1.8540 


1.8597 


1.8640 


1.8824 


1.8884 


1.8914 


1.8933 


u 


2.1876 


2.2100 


2.2213 


2.2289 


2.2346 


2.2571 


2.2642 


2.2677 


2.2698 



[5] and [7| the water amount is 4.05. The sensitivity on 7 is less than that of 
n, but is sufficient. Nevertheless, taking transient information into account, as 
given in the rows of Tables |6] and [7l will benefit the determination of 7 via this 
experimental set-up. 

5.7 Transient saturation determination 

In experiment 15.51 we have shown the low sensitivity of Mr on in an equi- 
librium analysis. One can assume this might improve in a transient analysis, as 
fiow around the interface will be determined by Ks- We investigate this now. 

To assure the sample is in part fully saturated, high rotational speeds are 
used. At fixed Ks we start from the equilibrium saturation at a; = 35 and let 
the centrifuge operate for lO^s (w 27.8hours) with lo = 40. After this time we 
change the rotational speed with the increment Aw = 5 and again operate for 
the same duration. This is continued up to rotational speed uj = 60, leading to 
5 end starting times j 10^, j = . . . , 4 and five end times j 10^, j = 1, . . . , 5. 

The same procedure is applied for a series of Ks values. Then we compare 
the values of Mr and Gc at the same time levels for different Ks- In Table [5] we 
present the values of Mr at time steps tj = 1000 + j IQ^ , j = 0, . . . , 4 (index j 
corresponds to the rows). The 9 columns correspond to the values Ks — kg 10~^ 
with ks = 0.24; 0.78; 1.32; 1.86; 2.4; 7.8; 13.2; 18.6 and 24. So, in this table 
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Figure 13: Head profiles corresponding with Fig. [T^ 



Table 9: Rotational momentum vs. kg under changing uj at the time level 





fcl 


fc2 


fc3 


/C4 


h 


fc'e 


fcy 


fc's 


fcg 


h 


0.9578 


0.9610 


0.9629 


0.9641 


0.9651 


0.9702 


0.9723 


0.9734 


0.9740 




1.2191 


1.2266 


1.2307 


1.2334 


1.2356 


1.2457 


1.2492 


1.2509 


1.2518 


h 


1.5127 


1.5253 


1.5319 


1.5364 


1.5398 


1.5542 


1.5587 


1.5609 


1.5622 


U 


1.8386 


1.8570 


1.8663 


1.8725 


1.8771 


1.8953 


1.9009 


1.9036 


1.9052 




2.1968 


2.2215 


2.2337 


2.2416 


2.2473 


2.2690 


2.2755 


2.2786 


2.2806 



Mr is computed 1000s after each change of w. The values of Mr at time level 
tj = j 10^, j = 1, . . . , 5 after each change of w is presented in Tabled 

In the Tables [TUl and [TTI we present the values of Gc which correspond to the 
values of Mr in Tables [8] and [9] (at the identical location) . 

Although the changes in the tables are clear and consistent, they remain low. 
Tablesl^land llll efFectivelv correspond to the equilibrium measurements, whereas 
the comparison with Tables [S] and [TU] gives an estimate of what is gained by a 
transient analysis. We can only conclude that for this set-up also a transient 
analysis to determine adds little improvements to the sensitivity. Note that 
for the parameters where an equilibrium analysis results in good sensitivity, 
a transient analysis is possible and can dramatically decrease the centrifuge 



Table 10: Center of gravity vs. k^ under changing u at the time level tj — 
1000 + j 10'^ 





fcl 


fc2 


fcs 




ks 


ke 


kj 


ks 


fcg 


to 


7.0485 


7.0532 


7.0580 


7.0591 


7.0608 


7.0745 


7.0826 


7.0885 


7.0946 


h 


7.0988 


7.1342 


7.1554 


7.1694 


7.1810 


7.2444 


7.2722 


7.2877 


7.2975 




7.1515 


7.2133 


7.2474 


7.2721 


7.2912 


7.3821 


7.4154 


7.4324 


7.4427 


h 


7.2011 


7.2840 


7.3288 


7.3600 


7.3836 


7.4866 


7.5206 


7.5376 


7.5480 




7.2462 


7.3460 


7.3981 


7.4335 


7.4596 


7.5660 


7.5993 


7.6159 


7.6261 
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Table 11: Center of gravity vs. fcs under changing lo at the time level tj = j 10"" 





fcl 


fc2 


fcs 


ki 


fcs 


fce 


kr 


fcs 


fcg 


h 


7.0965 


7.1279 


7.1469 


7.1591 


7.1693 


7.2235 


7.2462 


7.2579 


7.2644 




7.1496 


7.2085 


7.2419 


7.2644 


7.2823 


7.3672 


7.3971 


7.4115 


7.4197 


h 


7.1993 


7.2802 


7.3238 


7.3540 


7.3769 


7.4757 


7.5073 


7.5225 


7.5315 


U 


7.2447 


7.3429 


7.3941 


7.4288 


7.4544 


7.5578 


7.5894 


7.6048 


7.6140 


t5 


7.2859 


7.3977 


7.4539 


7.4910 


7.5178 


7.6210 


7.6518 


7.6671 


7.6763 
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Figure 14: Evolution of the water level in the reservoir, ks — 0.5; 1; 1.4; 2.4; 3.4 
and 4.4 

operating time. 

5.8 Saturated conductivity - 2 

In this experiment we demonstrate that the determination of model parameter 
Ks can be easily realized by means of very simple measurements obtained by 
centrifugation of a fully saturated sample with a water reservoir as presented 
in Section 31 In Fig. [U and (TS] the time evolution of the water level i{t) and 
of the rotational momentum Mr{t) is given. Note that in the computation of 
Mr{t) also a water collector placed at r G (rp + L, tq + L + d) is considered. 
The parameter values considered are Kg ks 10^^ with kg = 0.5; 1; 1.4; 2.4; 
3.4 and 4.4. The sensitivity of on ks will increase with higher values of u. 

In [5] £{t) for a fixed t G [0,Te] is used for the determination of Kg- From 
Fig. HH and US] we can deduce that also the values of Mr{t) (at a fixed t G [0, T^] 
) can be used for the determination of Ks- Adding transient data makes this 
determination still quicker. 

6 Acknowledgements* 

The first and the third author confirm financial support by the Slovak Research 
and Development Agency under contract APVV-0351-07. 



21 



2500 




1000' • • • • ' 1 

1000 2000 3000 4000 5000 6000 

X 



Figure 15: Evolution of the rotational momentum, ks — 0.5; 1; 1.4; 2.4; 3.4 and 
4.4 
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